Gravitational collapse of a Hagedorn fluid in Vaidya geometry 



T. Harko* 

Department of Physics, The University of Hong Kong, Pokfulam Road, Hong 

Kong, P. R. China. 

The gravitational collapse of a high- density null charged matter fluid, satisfying the Hagedorn 
equation of state, is considered in the framework of the Vaidya geometry. The general solution 
of the gravitational field equations can be obtained in an exact parametric form. The conditions 
for the formation of a naked singularity, as a result of the collapse of the compact object, are also 
investigated. For an appropriate choice of the arbitrary integration functions the null radial outgoing 
geodesic, originating from the shell focussing central singularity, admits one or more positive roots. 
Hence a collapsing Hagedorn fluid could end either as a black hole, or as a naked singularity. A 
possible astrophysical application of the model, to describe the energy source of gamma-ray bursts, 
is also considered. 
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CN . I. INTRODUCTION 

Investigating the final fate of the gravitational collapse of an initially regular distribution of matter, in the framework 
of the Einstein theory of gravitation, is one of the most active fields of research in contemporary general relativity. One 
^ I vifould like to know whether, and under what initial conditions, gravitational collapse results in black hole formation. 

One would also like to know if there are physical collapse solutions that lead to naked singularities. If found, such 
—i , solutions would be counterexamples of the cosmic censorship hypothesis, which states that curvature singularities in 
^ ■ asymptotically flat space-times are always shrouded by event horizons. 

I Roger Penrose [1] was the first to propose the idea, known as cosmic conjecture: does there exist a cosmic censor 
who forbids the occurrence of naked singularities, clothing each one in an absolute event horizon? This conjecture can 
be formulated in a strong sense (in a reasonable space-time we cannot have a naked singularity) or in a weak sense 
^ . (even if such singularities occur they are safely hidden behind an event horizon, and therefore cannot communicate 
■ with far-away observers). Since Penrose' s proposal, there have been various attempts to prove the conjecture (see [2] 
. and references therein). Unfortunately, no such attempts have been successful so far. 

Since, due to the complexity of the full Einstein equations, the general problem appears intractable, metrics with 
^H, special symmetries are used to construct gravitational collapse models. One such case is the two-dimensional reduction 
^ ' of general relativity obtained by imposing spherical symmetry. Even with this reduction, however, very few inhomo- 
geneous exact nonstatic solutions have been found. One well-known example is the Vaidya metric [3]. It describes the 
• • . gravitational field associated with the eikonal approximation of an isotropic fiow of unpolarized radiation, or, in other 
^ ' words, it represents a null fluid. It is asymptotically flat and it is employed in modeling the external field of radiating 
, stars and evaporating black holes. The second one is the Tolman-Bondi metric [4], which gives the gravitational field 
^ ' associated with dust matter and is frequently applied either in cosmological models or in describing the collapse of a 
5^ I star to a black-hole. Tolman-Bondi space-times embody the Schwarzschild solution, the Friedman universes and the 
Oppenheimer-Snyder collapse, as well as inhomogeneous expansions and collapses. 

At first sight these two metrics are completely different. Do the naked singularities that form in the collapse of 
null radiation and in the collapse of dust bear any relation with each other? Are there any features common to both 
solutions? And if this is the case, what are the implications for cosmic censorship? As shown by Lemos [5], the 
naked singularities which appear in Vaidya and Tolman-Bondi space-times are of the same nature. Various important 
features such as the degree of inhomogeneity of the collapse necessary to produce a naked singularity, the Cauchy 
horizon equation, the apparent horizon equation, the strength of the singularity and the stability of the space-time 
have a mutual correspondence in both metrics. For cosmic censorship, this result implies that if the shell-focusing 
singularities arising from the collapse of a null fluid are not artifacts of some (eikonal) approximation, then the shell- 
focusing singularities arising from the collapse of dust are also not artifacts (and vice versa) . Conversely, if the naked 
singularities are artifacts in one of them so are they in the other. 

Thus the Vaidya metric belongs to the Tolman-Bondi family. The most unbound case yields the Vaydia metric. 
Hence one expects that major features which might arise in one of the metrics will also appear in the other. One 
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example is the result that the strength in the Vaidya metric depends on the direction from which the geodesies enter 
the singularity [5]. 

Null fluids, are, in principle, easier to treat than matter fields. A null fluid is the eikonal approximation of a 
massless scalar field. Thus if one shows that the naked singularities arising in the Vaidya metric can be derived from 
more fundamental (massless) fields, then the naked singularities, which form in the Tolman-Bondi collapse, may also 
be derived from more fundamental (massive) fields. The same types of relations and conclusions hold for charged 
radiation and charged dust matter. The structure and properties of singularities in the gravitational collapse in Vaidya 
space-times have been analyzed, from different points of view, in [6]- [10]. 

Within the framework of various physical models, the spherical gravitational collapse has been analyzed in many 
papers. The role of initial density and velocity distributions towards determining the final outcome of spherical dust 
collapse and the causal structure of the singularity has been examined in terms of evolution of apparent horizon in 
[11], [12]. The collapse was described by the Tolman-Bondi metric with two free functions. The collapse can end in 
either a black hole or a naked singularity. The occurrence and nature of naked singularities in the Szekeres space-times 
have been investigated in [13]. These space-times represent irrotational dust. They do not have any Killing vectors 
and they are the generalizations of the Tolman-Bondi space-times. There also exist naked singularities that satisfy 
both the limiting focusing condition and the strong limiting focusing conditions. The relevance of the initial state 
of a collapsing dust cloud towards determining its final fate in the course of a continuing gravitational collapse has 
been considered in [14]. Given any arbitrary matter distribution for the cloud in the initial epoch, there is always the 
freedom to chose the rest of the initial data, namely the initial velocities of the collapsing spherical shells, so that, 
depending on this choice, the collapse could result either in a black hole or a naked singularity. Thus, given the initial 
density profile, to achieve the desired end state of the gravitational collapse one has to give a suitable initial velocity 
of the cloud. The expression for the expansion of outgoing null geodesies in spherical dust collapse has been derived in 
[15] The limiting values of the expansion in the approach to singularity formation have been computed. Using these 
results one can show that the covered, as well as the naked singularity solutions arising in spherical dust collapse, are 
stable under small changes in the equation of state. 

The growth of the Weyl curvature is examined in two examples of naked singularity formation in spherical gravi- 
tational collapse-dust and Vaidya space-time- in [16]. The Weyl scalar diverges along outgoing radial null geodesies 
as they meet the naked singularity in the past. Although general relativity admits naked singularities arising from 
gravitational collapse, the second law of thermodynamics could forbid their occurrence in nature. A simple model for 
a corona of a neutrino-radiating star showing critical behavior is presented in [17]. The conditions for the existence or 
absence of a bounce (explosion) are discussed. The charged Vaidya metric was extended to cover all of the space-time 
in [18] and the Penrose diagram for the formation and evaporation of a charged black hole obtained. The covariant 
equations characterizing the strength of a singularity in spherical symmetry and a slight modification to the definition 
of singularity strength have been derived in [19]. The idea of probing naked space-times singularities with waves 
rather than with particles has been proposed in [20]. For some space-times the classical singularity becomes regular 
if probed with waves, while stronger classical singularities remain singular. 

In order to obtain the energy-momentum tensor for the collapse of a null fluid an inverted approach was proposed 
by Husain [21]. First the stress-energy momentum tensor is determined from the metric. Then the equation of state 
and the dominant energy condition are imposed on its cigc;nvalues. This leads to an equation for the metric function. 
The precise form of the stress-energy tensor is then displayed. By using this approach several classes of solutions 
describing the collapse of a null fluid, satisfying barotropic and polytropic type equations of state, have been obtained. 
In the framework of the same approach a large class of solutions, including Type II fluids, and which includes most 
of the known solutions of the Einstein field equations, has been derived, in four dimensions by Wang and Wu [22], 
and in A'' > 4 dimensions by Villas da Rocha [23]. The Vaidya radiating metric has been extended to include both a 
radiation field and a string fluid by Glass and Krisch [24], [25] and by Govinder and Govender [26]. 

When nuclear matter is squeezed to a sufficiently high density, a phase transition takes place and neutron matter 
converts into three-flavor (strange) quark matter, which is due to the fact that strange matter may be more stable than 
nuclear matter. The collapse of the quark fluid, described by the bag model equation of state p = {p — AB) /3, with 
B =constant, has been studied by Harko and Cheng [27] and the conditions for the formation of a naked singularity 
have been obtained. The obtained solution has been generalized to arbitrary space-time dimensions and to a more 
general linear equation of state by Ghosh and Dadhich [28] , [29] . 

In 1965 Hagedorn [30] postulated that for large masses m the spectrum of hadrons p (m) grows exponentially, 
p{m) ^ exp(m/Tff), where Th, the Hagedorn temperature, is a scale parameter. The hypothesis was based on the 
observation that at some point a further increase of energy in proton-proton and proton-antiproton collisions no longer 
raises the temperature of the formed fireball, but results in more and more particles being produced. Thus there is a 
maximum temperature Th that a hadronic system can achieve. The statistical model of the hadrons has been used 
to obtain a description of dense matter at densities exceeding nuclear density. The Hagedorn phase also arises in 
theories containing fundamental strings, because they have a large number of internal degrees of freedom [31]. As a 
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result of the existence of many oscillator modes the density of states grows exponentially with single string energy. 
Thermodynamical quantities, such as the entropy, diverge at the Hagedorn temperature. If one considers an ensemble 
of weekly interacting strings at finite temperature, this behavior of the density of states is thought to lead either to a 
limiting temperature or a phase transition, in which the string configuration changes to one which is dominated by a 
single long string [32] . The high density Hagedorn phase of matter has been extensively used in cosmology to describe 
the very early phases of the evolution of the Universe [33]- [36]. 

It is the purpose of the present paper to study the spherically symmetric gravitational collapse of the charged 
matter in the Hagedorn phase. In order to simplify the mathematical formalism we adopt the assumption that the 
high density fluid moves along the null geodesies of a Vaidya type space-time. The Vaidya geometry, also permitting 
the incorporation of the effects of the radiation, offers a more realistic background than static geometries, where all 
back reaction is ignored. By adopting the Hagedorn equation of state for dense matter, the general solution of the 
field equations can be obtained in an exact form. For the sake of comparison we also consider the collapse of matter 
described by the stiff Zeldovich equation of state. The conditions of formation of naked singularities are obtained in 
both cases. 

The present paper is organized as follows. The field equations describing the collapse of a Hagedorn fluid are 
written down in Section II. The general solution of the field equations is presented in Section HI. The equation of the 
null outgoing geodesies and the conditions of the formation of the naked singularities are discussed in Section IV. An 
astrophysical application of the formalism to explain the gamma ray bursts energy emissions is described in Section 
V. In Section VI we discuss and conclude our results. 



II. GEOMETRY AND FIELD EQUATIONS 



In ingoing Bondi coordinates {u, r, 9, ip) and with advanced Eddington time coordinate u = t + r (with r > the 
radial coordinate and r decreasing towards the future) the line element describing the radial collapse of a coherent 
stream of matter can be represented in the form [21], [27] 



2m {u, r) 



dvT + 2dudr + 



+ sin^ 9dip^ 



(1) 



m {u, r) is the mass function and gives the gravitational mass within a given radius r. In the following we use the 
natural system of units with SttG = c = 1. 

The matter energy-momentum tensor can be written in the form [21], [22] 



where 



is the component of the matter field that moves along the null hypcrsurfaces u = const., 

= iP + P) (If^n^ + I'^n^) +pg^u, 
represents the energy-momentum tensor of the collapsing matter and 



(2) 

(3) 
(4) 



(5) 



is the electromagnetic contribution. and are two null- vectors given by Ifj, = and 



2m(u,r) 



S^y,. so that Z„P 

(m) ■ " 



an" = and ^^n" = — 1 (with S^^'j^^ the Kronecker symbol) [21], [22]. The energy density and 
pressure in Eq. (4) have been obtained by diagonalizing the energy- momentum tensor obtained from the metric [21]. 
The electromagnetic tensor F^^, obeys the Maxwell equations [37] 



dxt" 



(6) 



-g dxi^ 



(7) 



3 



Without any loss of generality the electromagnetic vector potential can be chosen as [6], [7] 

with q{u) being an arbitrary integration function. From the Maxwell Eqs. (6) and (7), it follows that the only 
non-vanishing components of F^;/ are Fru = ~Fur = ^^r- and, consequently, 

E''^ = ^diag {-1,1, -1,1). (9) 

For the energy-momentum tensor (2) the gravitational field equations take the form [27] 

1 dm (u, r) 1 , , „, 

+ (11) 

-- Q,2 ' =p{u,r)+^-y. (12) 
The stress-energy tensor (4) satisfies the dominant energy condition if the following three conditions are met: 

p>0,p>p,Tabw''w^>0, (13) 

where w"' is an arbitrary time-like four-vector. The first two of these conditions imply that ^ > and ^^J- < 0. 
The former just says that the mass function either increases with r or is a constant, which is a natural physical 
requirement on it. To satisfy the first two of the dominant energy conditions one must impose an equation of state 
for the collapsing matter. 

Usually two different equations of state are used for the description of matter at extremely high densities. One of 
the most widely investigated case is the so-called causal limit of the linear barotropic equation of state p = (7 — 1) p, 
7 =constant, corresponding to 7 = 2, or the Zeldovich stiff fluid equation of state p = p . 

The Zeldovich equation of state, valid for densities significantly higher than nuclear densities, p > lOpnuc, with 
Pnuc = 10^^ g/cm^, can be obtained by constructing a relativistic Lagrangian that allows bare nucleons to interact 
attractively via scalar meson exchange and repulsively via the exchange of a more massive vector meson [38] . In the 
non-relativistic limit both the quantum and classical theories yield Yukawa-type potentials. At the highest densities 
the vector meson exchange dominates and by using a mean field approximation one can show that in the extreme limit 
of infinite densities the pressure tends to the energy density, p ^ p. In this limit the sound speed Cs = y^dp/dp 1, 
and hence this equation of state satisfies the causality condition, with the speed of sound less than the speed of light 

An alternative approach to the equation of state at ultra- high densities is based on the assumption that a whole host 
of baryonic resonant states arise at high densities. In the Hagedorn model the baryon resonance mass spectrum, that 
is the increase in the number of species of particles with mass between m and m + dm, is given by dN = N{m)dm ~ 

oexp {m/Tn) / {m^ + m^)' dm, where N{m)dm is the number of resonances between mass m and m + dm. Fitting to 
the existent experimental data on baryon resonances show that t7?.o = 500 MeV and a = 2.63 x 10'* MeV'^/^ [39]. If dN 
increased any faster as m ^ 00 than in the above formula, the partition function would not converge. Also the partition 
function converges only if the temperature of the system is less than Th- Thus Th, the Hagedorn temperature, is the 
effective highest temperature for any system. By interpreting the transverse momentum distribution of secondaries in 
very high energy collisions in terms of the model one obtains Th ^ 150—190 MeV and Z = 5/4 [39]. The corresponding 
equation of state for the matter is 

P = Po + /9oln— , (14) 

Po 

where po = 0.314 x lO^'^g/cuv^ and po = 1.253 x 10*^ g/cm"^ [40], [41]. The velocity of sound in this type of matter 
is Cs = yPoJp ■ For the Hagedorn equation of state the speed of sound has the property Cg ^ for p/po 00, in 
striking contrast with the mean field theory approach, leading to the Zeldovich equation of state, in which Cg — > 1. 
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The Hagcdorn equation of state creates new " particles" continuously with increasing density, rather than enlarging 
the Fermi sea of a single species. The equation of state (14) could be valid asymptotically for densities greater than 
about 10 times the nuclear density p„ = 2 x 10^^ g/cm^. The vacuum boundary of the initial matter distribution is 
defined by the equation p = 0, condition corresponding to a surface density ps = po&~^^^ ~ 0.778 x 10^'' g/cm^, two 
times smaller in magnitude as the nuclear density. This condition also defines the physical radius of the initial matter 
distribution and defines a boundary for the null fluid. 

Hence, as a possible physical model to describe high-density matter in the final stages of the gravitational collapse 
we shall adopt the Hagedorn equation of state. In the following Section we present the general solution of the 
gravitational field equations for the null Hagedorn fluid. 



III. SPHERICAL COLLAPSE OF THE HAGEDORN NULL FLUID 



From Eq. (11) we immediately obtain 

Po 



^dm{u,r) q^{u) 
dr 



(15) 



With the use of the Hagedorn equation of state Eq. (14) and of Eq. (15), Eq. (12) can be written in the form 



d'^m{u,r) q^iu) 



Qj.2 



= Po+ Po In 



por^ 



^dm{u,r) q^{u) 
dr 



By introducing a new variable 



w (u, r) = In — = In 
Po 



por^ 



dm{u,r) q^{u) 



dr 



it is easy to show that 



1 
Por 



dr^ r^ 



1 dw 



Therefore Eq. (16) becomes 



dw „ a + w + e 

r— — = -2 

dr 

where a = po/po ^ 0.25. By integrating Eq. (19) we obtain 



C{u) 



= exp 



2/a 



+ w + e'" 







j = exp 





(16) 



(17) 



(18) 



(19) 



(20) 



where C{u) is an arbitrary integration function and we denoted F{w) = J dw. Eq. (20) formally defines w 

as a function of r] = r/C{u), w = H {rf). The variation of the mass function is described by the equation 



2^^ = ^ + por^e-, 
dr 



having the general solution given by 



2m(n,r) = D{u) - ^e^-(-) - («) / (^^ " I^N) 
^ ' ' ^ ' C{u) 2 ^ 'J a + w + e"" 

where D (u) is an arbitrary integration function. In the following we shall also denote 



dw, 



K{w) 



r exp (2w - ^F{w)) 
J a- 



dw. 



' + w + e'^ 

The density and the pressure of the collapsing null Hagedorn fluid are given by 



(21) 



(22) 



(23) 
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p (w) = po exp (w) ,p{w) = po + pow. (24) 
The energy density /x of the radiation moving along the u =constant null hypersurfaces is given by 

'^'^ ' ' C^{u) C^{u) 2 ^ ^ ^ ^ 2 ^ ^a + w + e'^ dr] ^ ^ 

where a dot denotes the derivative with respect to u. 

The function r/C{u) can be represented as a power series of the parameter w in the form 

r 2w Up' nfl liw'^ 367w^ ^r, ^ . ^ 

1-^ + ^- ^ + ^- -T^ + 0'[w]. (26) 



C{u) 5 5 5 60 1500 

while the function K{w) has the power series representation 

, 4w 8u;2 52^3 59u;* 351u;^ , 

K{w)k \ 1 hO M- (27) 

^ ^ 5 25 25 125 625 ^ ^ ^ ' 

In the limit of large densities p — > oo, which also implies w — > oo, Eq. (19) becomes 



with the general solution given by 



M;(w,r) «ln^^. (29) 



Therefore in the asymptotic limit of very high densities, corresponding to r ^ 0, the gravitational collapse of the 
Hagedorn null fluid is described by the equations 

2m{u,r)KD{u)+pQC{u)r-^^, (30) 

, , Ciu) , , , C[u) 

p{u,r)'» po—^,p{u,r)KpQ + po\n—^, (31) 

/ N D C 2qq 
/.(t*,r)«^+Po--^. (32) 

In order to find the behavior of the solution in the opposite limit of large, but finite r, wc note first that the 
boundary of the Hagedorn type matter distribution is defined by the equation p= 0, corresponding to a value of the 
parameter w = Wg = —po/po ~ —1/4. Values of w < Wg lead to the unphysical situation of negative pressure matter 
distribution, p < for w < Ws. In order to find the behavior of the solution near Wg, wc represent w{u,r) in the 
form w{u,r) = Ws + wi{u,r), with wi {u,r) a small perturbation satisfying the condition \wi {u,r)\ << \ws\, ^u,r. 
Substituting into Eq. (19) gives 

^ dr ^ ei«.e«'i ~ s(l + u;i) ' ^ ' 



where we denoted s = exp (ws) « exp(— 1/4) w 0.778. Hence for large r the solution of Eq. (33) is given by 

= -21nr + C"(w), (34) 



s sins s , 
-wi + -o + -o In 



1 + s {l + sY (1 + s) 



1 H wi 



where C'{u) is an arbitrary integration function. Since we have assumed that wi {u, r) is small, we can neglect in Eq. 
(34) the term containing the logarithmic function. Hence near the vacuum boundary of the high density Hagedorn 
fluid distribution we obtain 
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1 C(w) 
w K. too +ln^^, 



(35) 



where wq = Wss/{1 + s) —0.10, fc = (1 + s) /s = 2.284 and C{u) is an arbitrary, u-dependent integration function. 
Therefore in the limit r — > oo the general solution of the gravitational field equations can be approximated by 



2m {u, r) « D{u) — 



3- 2A: 



C{u) 



3-2fc 



(36) 



p(.,r)«,oe-^«.0.9po^. 



/ ^ , Ciu) / , C{u) 

p [u, r)^po + PqWq + po In w po I 0.15 + In 



(37) 
(38) 



Mn,.)«2^-?^ + ^C(.)/--. 
b{u) 2q{u)q{u) ^.^^C{u) 

2-^:2 ;:3 0-58po^3:56, 

where D(u) is an arbitrary integration function. 

The variation of the ratio r/C{u) as a function of w is represented, by using Eq. (20), in Fig. 1. 



(39) 



u 



1 

0.8 
0.6 
0.4 
0.2 



01234567 
w 

FIG. 1. Variation of the ratio r/C{u) as a function of the parameter w. 

As one can see from the figure, in the limit w ^ Wg = —0.25, the ratio r/C{u) tends to 1, while for large w, w ^ og, 
that is, in the limit of very high densities, r/C{u) tends to zero. Values oi w < Ws = —0.25 lead to the violation of 
the dominant energy condition p > 0. The function K{w) is represented, as a function of r/C{u), in Fig. 2. 
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5 
1^ 




0.4 0.6 

r/C(u) 

FIG. 2. Variation of K{w) as a function of r/C{u). 



For the sake of comparison wc shall also present the solution of the gravitational field equations corresponding to 
the collapse of the charged null Zeldovich fluid with p = p. In this case Eqs. (11) and (12) give immediately 



with the general solution given by 



2 dm {u, r) 1 d'^m {u, r) 



dr 



m{u,r) = B{u) 



(40) 



(41) 



where A[u) and B(u) are arbitrary integration functions. 

The other physical quantities characterizing the collapsing Zeldovich fluid are given by 



p{u,r) =p{u,r) 



2^2 (u) - g2(y) 



(42) 



, , 2 fdB 2 ^, .dA 



du 



du 



(43) 



In all these cases the electromagnetic current follows from the Maxwell Eq. (7) and can be generally represented as 



1 dq{u) 
Airr'^ du 



(44) 



The energy-momentum tensor of the mixture of fluids under consideration belongs to the Type II fluids [21]. The 

energy conditions arc the weak, strong and dominant energy conditions fJ.>0,p>0,p>0,p>p>0 and can be 
satisfied for both Hagedorn and Zeldovich fluids by appropriately choosing the arbitrary functions A{u), B{u), C{u) 
and D{u) that characterize the injection and initial distribution of the mass and q{u) that describes the variation 
of the charge. For the Hagedorn fluid the conditions p > and p > urc trivially satisfied for the values of the 
parameter w in the range w € (ws,oo), corresponding to the range of densities of p € (ps,oo), with Ws = —0.25 and 
Ps = 0.78 X 10^** g/cm^. Due to the choice of the parameters in the equation of state the dominant energy condition 
p> p is also satisfied. The condition > is equivalent to ^ > 0, and in the approximation of high densities leads 
to 



dPju) dC{u) ^ dq{u) 1 
du du ^ ~ ^ ^ du r' 



(45) 



imposing a simultaneous constraint on all thrc;e functions C{u), D{u) and q{u). For small values of r and for a charged 
collapsing high density fiuid, the right hand side of Eq. (45) dominates and this energy condition could not hold. 



One possibility to satisfy Eq. (45) for all r is to assume that the function q(u) behaves so that 



dq^(u) 
du 



for r ^ 0. 



This means that the charge in the singular point r = is constant for all times. Alternatively, we may suppose that 
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at extremely small radii matter is generated so as to satisfy the energy eondition. For neutral q = matter Eq. (45) 

> and ^ 

du du 



is easily satisfied by choosing i^^M > and ^^^IhI > q jn the case of the Zeldovich fluid the energy conditions are 



satisfied by choosing the functions A, B and q so that A'^{u) > q^{u)/2 and B{u) > 2A{v)A [u) jr. 

The radii of the apparent horizon of the metric (1) are given by the solution of the equation 2m = r. If 
lim„^ooC'(M) = Co = const., limi,^oo -D(u) = -Do = const, and lim„^oo = go = const., then the algebraic 
equation determining the radii of the apparent horizons in the case of the Hagedorn fluid is 

Coe-i^(-) + ^CiK{w) + ^ei^M = Do, (46) 
^ Co 

which in general may have multiple solutions. 

In the case of the collapse of the Zeldovich fluid the radii of the apparent horizon are given by the solutions of the 
equation 

2Bo -2^=r, (47) 

r 

where Aq = limi,^oo C(m) =constant and Bq = lim„^oo D{u) =constant. The radii of the apparent horizon are 

ri,2 = Bo± ^JbI - 2Al (48) 

The singularities of the matter filled Vaidya space-time can be recognized from the behavior of the energy density 
and curvature scalars like e.g. RapR"^ and Ra^-ysR"''^^^ , given in the high density limit for the collapsing Hagedorn 
fluid by 

i?a/3-R"^ = ^ [rVoC'(^) + 2r''poC{u)q\u) + 2q\u)] , (49) 

= ^ r-VgC^(u) - 2r^poC{u)D{u) + Zr-^p-^ju) - 2r^poCiu)q\u) + 12rD{u)qHu) + lAq^u) ^ ^^^^ 
and which diverge for r ^ 0. 

IV. OUTGOING RADIAL NULL GEODESICS EQUATION 

The central shell-fociissing singularity (i.e. that occurring at r = 0) is naked if the radial null-geodesic equation 
admits one or more positive real roots Xq [2]. In the case of the pure Vaidya space-time it has been shown that for 
a linear mass function 2m{u) = Xu the singularity at r = 0, tt = is naked for A < | [42]. Hence it is important 
to investigate whether the gravitational collapse of high density matter described by the Hagedorn and Zeldovich 
equations of state could result in the formation of naked singularities. 

We consider first the case of the gravitational collapse of the Hagedorn fluid. In order to simplify the calculations 
we chose some simple particular expressions for the functions C{u), Diu) and q{u), e.g. C{ti) = agu, D{u) = P^yu 
and q{u) = qou, with ao > 0, /3o > and qo > constants. With this choice the equation of the radially outgoing, 
future- directed null geodesic originating at the singularity can be written as 

^ = ? (51) 

l-(3o{f)+qUT)' + ^^Hf)KM 

For the geodesic tangent to be uniquely deflned and to exist at the singular point r = 0,u = of Eq. (51) the 
following condition must hold [2] : 

lim -= lim ^=Xo. (52) 
u,r—^0 r M,r— >o dr 

When the limit exists and Xq is real and positive, there is a future directed, non-space-like geodesic originating 
from r = 0, M = 0. In this case the singularity will be, at least, locally naked. 
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Since the function K[w) is finite for all w = H {r/C{u)), and wc assume that Imiu.r-^ou/r is also finite, it follows 
that limu^r^ou^ (7) K{w) = 0. Therefore it follows that for the null geodesic Eq. (51) condition (52) leads to the 
following third order algebraic equation: 

(Zo'^o - /3o^o +Xo-2 = 0. (53) 

In the case of a neutral Hagcdorn fluid with q(u) = 0, Eq. (53) reduces to a second order algebraic equation with 
two roots, Xoi,2 = (1 ± y/1 — 8/3o) /2/3o- Therefore the condition for the formation of a naked singularity is /3o < 1/8. 
For Pq > 1/8, as a result of the collapse of the Hagcdorn type fluid, a black hole is formed. 

For q(u) ^ the condition of the existence of at least one real solution of Eq. (53) is 

lOSq^o + 4^0 + 8/3o' > Po + 36/3oao • (54) 

Therefore, by appropriately choosing the constants go and /Jo, it is always possible to construct a positive solution 
ofEq. (53). 

In the case of the collapse of the Zeldovicli fluid and by assuming for the arbitrary u-dependent integration functions 
the form A{u) = a^u and B{u) = bou, with ao, bo non-negative constants, the equation of the radially outgoing, future- 
directed null geodesic originating at the singularity is 

— = ^ (55) 

dr i_6„(H)+«2(H)- 

The algebraic condition for the formation of a naked singularity is given by the requirement that the equation 

alX^ - 60^0 +^0-2 = 0, (56) 

has at least one positive root Xq > 0. The condition that the above equation has at least one real root can be written 

as 

a| ^ l-8bo 

bl - 108a2-366o + 4' ^ ' 

This condition can be satisfied, for example, by choosing 60 < 1/8 and ao > l/-\/216. Therefore, as in the case of 
the charged Hagedorn fluid, it is always possible to construct a positive solution of Eq. (56). 



V. COLLAPSING HAGEDORN MATTER-A POSSIBLE SOURCE OF GAMMA-RAY BURSTS 

Gamma-ray bursts (GRBs) are cosmic gamma ray emissions with typical fluxes of the order of 10"^ to 5 x 10~^erg 
cm~^ with the rise time as low as 10~^ s and the duration of bursts from 10~^ to 10^ s [43]. The distribution of 
the bursts is isotropic and they are believed to have a cosmological origin, recent observations suggesting that GRBs 
might originate at extra-galactic distances [43] . The large inferred distances imply isotropic energy losses as large as 
3 X 10^3 gj.g for GRB 971214 and 3.4 x lO^" erg for GRB 990123 [44]. In the present Section we shall use CGS units. 

The widely accepted interpretation of the phenomenology of 7-ray bursts is that the observable effects are due to 
the dissipation of the kinetic energy of a relativistically expanding fireball, whose primal cause is not yet known [43]. 

The proposed models for the energy source involve merger of binary neutron stars [45], capture of neutron stars 
by black holes [46], differentially rotating neutron stars [47] or neutron star-quark star conversion [48] etc. However, 
the most popular model involves the violent formation of an approximately one solar mass black hole, surrounded 
by a similarly massive debris torus. The formation of the black hole and debris torus may take place through the 
coalescence of a compact binary or the collapse of a quickly rotating massive stellar core [49] . There are still many 
open problems concerning GRBs, from which the most important is the problem of the source of the large energy 
emission during the bursts. On the other hand naked singularities as sources of 7-ray bursts have also been proposed 
[50], [51], [52]. The fact that explosive radiation can be emitted during the gravitational collapse to a naked singularity 
of a dust ball has also been pointed out recently by [53] . 

As an astrophysical application of the Hagedorn fluid collapse in the Vaidya geometry we consider the possibility 
that gamma - ray bursts could in fact be energy emission during the collapse of a neutral, 9 = 0, high density matter 
in the Hagedorn phase, ending with the formation of a naked singularity. An estimation of the energy emitted during 
the collapse shows that it is of the same order of magnitude as the one measured during 7— ray bursts. Hence this 
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mechanism could provide a valuable explanation for this phenomenon, also opening the possibility of the observational 
investigation of the astrophysical properties of the naked singularities. 

The arbitrary integration functions C{u) and D{u) appearing in Eq. (22), describe the injection and initial distri- 
bution of the mass, respectively. Their exact mathematical form cannot be obtained in the framework of the present 
formalism. The only requirements that the functional form of these functions must satisfy are the dominant energy 
conditions and the condition that the equation of the radially outgoing, future-directed null geodesies, originating 
at the singularity r = 0, be uniquely defined and exist at the singular point r = and has at least one positive 
root. A positive root of the geodesic equation leads to the possibility of formation of a naked singularity as a result 
of the collapse. Some simple integration functions satisfying the condition on the geodesic equation, as well as the 
dominant energy conditions, are, for example, C(u) = ao^, D{u) = Pqu, with ao > 0, /3o > constants. With this 
choice a naked singularity may form during the collapse of the Hagedorn fluid if and only if the algebraic equation 
Po^o — + 2 = has at least one positive root Xq > 0. This condition requires /3o < 1/8. Of course many other 
choices are also possible. 

Naked singularities are gravitational singularities that are not covered by a horizon. Near the singularity the space- 
time curvature and the gravitational tidal forces grow very strongly. During the collapse naked singularities could 
emit powerful bursts of radiation visible to a distant external observer situated far away from the sight of the collapse 
[50]. 

As a first physical parameter we need to estimate is the time tff necessary for a matter element at the surface of 
the star to reach the center r = of the collapsed object, as a function of the mass at the center. This can be done 
by evaluating the time derivative of the mass, given by Eq. (30) with q = 0, at the center of the naked singularity: 



dm / dD{u 



dt 2G V du 

By integrating the above equation from t = to tff we obtain 



(58) 

r=0 



f'" fdDiu) , 

M 1^=0= ^ / — dt, (59) 



2G In \ du 



r=0 



where we assumed rn(0) |r=o= and denoted m{tff) \r=o= M |r-=o= 0. In order to find an explicit expression for 
tff we shall approximate the integral in Eq. (59) by using the first mean value theorem for integrals, which states 

that for any arbitrary function f(t), f{t)dt = {b — a)f{c), where a < c < b. f{c) is called the average value of /. 
Hence, in the following we shall approximate the derivatives of the arbitrary the arbitrary integration functions and 
the functions themselves by their average values. Therefore from Eq. (59) we obtain 

2G 

tff = -^M \r=o, (60) 
where we denoted by /3o > the average value of the function {dD{u)/du)^^Q, f3o = [dD{t)/d{t)]^^g, where < 6 < 

teff- 

The initial mass distribution of the Hagedorn fluid can be obtained from the study of the mass distribution at the 
initial moment t = 0: 



dm c? 
rfT '*=°= 2G 



dD(u) , dC(u) , 

du du 



(61) 



By integrating this equation, and using again the flrst mean value theorem, it follows that the initial mass profile 
of the Hagedorn fluid can be represented as 

m(r) |t=o= ^aor, (62) 

where ao ^ is the average value of the r-dependent function h{r) — [dD{u) / du + por{dC{u) / du)] \t=Q, with the 
function h{r) estimated at some point r = a,0<(j<R, with R the radius of the Hagedorn fluid distribution. Hence 
ao = h{a). From this equation we can see that our approach implies a linear profile of the initial mass distribution 
of the star, ao, the average value of the spatial derivatives of the integration functions D{u) and C{u) at f = 0, is 
completely determined by the initial conditions, that is by the total initial mass of the star M \t=o and by its initial 
radius R via the relation 
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ao = (63) 

In order to calculate the energy emitted during the Hagedorn fluid collapse we shall admit that the luminosity of 
the collapsing object should not exceed the rate of collapsing matter energy [52]. The variation of the mass of the 
high density fluid during collapse is given by 



^ 2Gm{u,r) ^ ^^^^^^ ^ porACiu) + poCiu)Ar] 



dD{u) ^ dCiu) ^ ^ ^ ■ 

du du 



(64) 



The variation of the advanced time coordinate u is given by the approximate expression Au = cAt + Ar = 
cAt (l + ~ cAt (l + where we defined Vf = ^ as being the speed of the collapsing fluid as measured by 
a local observer. Then we obtain 



^ 2Gm{u,r) 



du du c \ c 



: (l + ^) At. (65) 



We assume that the energy emission occurs mainly from a small region near the center of the naked singularity. 
Hence we shall evaluate Eq. (65) near r very close to zero. Then from Eq. (58) we can roughly approximate 
dD{u)/du K, (2G/C-'') [dm/dt) k, 2G/c^ (M/tf j) = 13^. Since for a small r, C {ct + r) w C(ct), wc can also neglect 
in this limit the term p^r {dC (u) / du) . From the field Eq. (11) and from Eq. (24) it immediately follows that 
2{dm{u,r) /dr) = (SttG / c^) r'^ p « poC{u). Taking into account Eq. (62) we immediately find poC{u) w Anao ~ 
(SttGM/c^R^. The velocity Vf of the collapsing matter can be obtained from vj » R/t/j ~ {c?RPo)/ {2GM |r=o)- 

With the use of the previous results we obtain for the total energy emitted per unit time during the collapse of a 
Hagedorn fluid to a naked singularity the expression: 



AEr 1 



At 2 



/3o+47rao^(l + ^)"'j (l + ^) ^ erg s"!. (66) 



Hence we have expressed the total energy of the radiation, which could be emitted from a naked singularity, in 
terms of the average values of the derivatives of the arbitrary integration functions. By assuming an emission time of 
the order of At ^ 10~^ s and a velocity of the collapsing matter of the order Vf = & x lO^cm/s we obtain a value of 

the emitted energy AE^ sa 1.85 x /3o + Anao^ (l + ■^) ^ x lO^^erg. Therefore the energy which could be emitted 

during the collapse of the Hagedorn fluid to a naked singularity could exceed in magnitude the energy AEj « 4 x 10^^ 
erg, observed for GRB990123 [44]. Of course the exact value of the energy depends on the exact numerical values of 
the constants ao and (3o- By comparing Eq. (66) with the observational data it is possible to determine the values of 

ao and Po- 

Under the effect of the collapse the compact object will heat up to a temperature of the order of T ~ Th ~ 10^'^K, 
higher than that occurring in super nova explosions. Since Th cannot be exceeded, the gravitational energy of the 
compact object is converted into new particles, which will be generated during the collapse. Most of the newly created 
particle will decay via weak interaction processes. As a result a neutrinosphere will form around the naked singularity. 
Therefore the main energy loss mechanism of the super-heated collapsed astrophysical object would be neutrino 
radiation. The neutrinos and anti neutrinos interact with protons and neutrons via the URCA processes n + ^ 

p+e~ , p+Ve n + e+. At temperatures higher than the nuclear Fermi temperature kT^'^ — ( j (2^) (y) > 

that can also be expressed in the form Tf{^ = 1.47 x 10~^pg^^ [49], the integrated optical neutrino depth is unity. 
Hence the deposition energy can be estimated as E {1 — e~'^) AE^ ergs. The process 7 + 7 ^ e+ + e~ will generate 
a fireball that will expand outward. The expanding shell interacts with the inter stellar medium surrounding the 
Hagedorn type naked singularity, and the kinetic energy is finally radiated through non-thermal processes in shocks 
[48]. 

The energy released during the collapse of the Hagedorn and radiation fluids into a naked singularity has the same 
order of magnitude as that observed in the case of gamma ray bursts. This strongly suggests the possibility that 
gamma ray bursts could be massive compact objects, formed from a Hagedorn fluid, collapsing to a naked singularity 
in a cosmological environment. 

VI. CONCLUSIONS 

In the present paper we have considered the collapse of a Hagedorn type fluid in the Vaidya geometry. 
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The exact solution obtained represents the generalization, for Hagcdorn type matter, of the collapsing solutions 
previously obtained by Vaidya [3], Bonnor and Vaidya [6], Lake and Zannias [7] and Husain [21]. From a mathematical 
point of view the solution is represented in a parametric form. It satisfies all the energy conditions and consequently 
describes the collapse to a singular state. The possible occurrence of a central naked singularity has also been 
investigated and it has been shown that, at least for a particular choice of the parameters, a naked singularity is 
formed. Depending on the initial distribution of density and velocity and on the constitutive nature of the collapsing 
matter, either a black hole or a naked singularity is formed. The values of the parameters in the solution (22) and 
(24) determine which of these possibilities occurs. The solution describing the collapse of both the Hagedorn and 
Zeldovich matter is asymptotically flat, but this condition does not play any significant role in the formation of the 
naked singularity. 

The Hagedorn type matter may reside as a permanent component of neutron stars core at high temperatures or 

densities and form stable compact stellar objects. In fact, from a physical point of view, it seems that the high density 
limit for the equation of state described by the Hagedorn equation of state is one of the best and more realistic 
candidates for the study of properties of collapsing objects. It also serves to illustrate the much richer interplay that 
can occur among particle physics and general relativity, when more involved quantum field theoretical models are 
considered. 

As a possible simple astrophysical application of this collapsing solution we have considered the possibility that 
gamma ray bursts could be energy emission during the collapse of a high density star, ending in the formation of a 
naked singularity. The radiated energy during this process could be as high as lO^^ergs. Thus the naked singularity 
explosion could be a candidate for the central engine of a gamma ray burst. 
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